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Abstract 

This paper examines an averaging technique in which the nonlinear 
flux term is expanded and the convective velocities are passed through 
a low-pass filter. It is the intent that this modification to the nonlinear 
flux terms will result in an inviscid regularization of the homentropic Euler 
equations and Euler equations. The physical motivation for this technique 
is presented and a general method is derived, which is then applied to 
the homentropic Euler equations and Euler equations. These modified 
equations arc then examined, discovering that they share the conservative 
properties and traveling wave solutions with the original equations. As 
the averaging is diminished it is proven that the solutions converge to 
weak solutions of the original equations. Finally, numerical simulations 
are conducted finding that the regularized equations appear smooth and 
capture the relevant behavior in shock tube simulations. 

1 Introduction 

This paper examines a modification of the homentropic Euler and Euler equa- 
tions where flux term is expanded and the convective velocities are passed 
through a low-pass filter. It is the intent that this modification to the non- 
linear flux terms will result in an inviscid regularization of the homentropic 
Euler equations and Euler equations. The ultimate goal with this regulariza- 
tion is to develop a proper modeling of small scale behavior so that both shocks 
and turbulence can be captured in one comprehensive technique. This paper 
primarily addresses the shock capturing capabilities of the developed technique. 

The inspiration for this new technique comes from recent work done on a 
similarly modified Burgers equation. The Burgers equation, Ut + uUx = 0, is 
considered a simplistic model of compressible flow which forms shocks readily. 
Classically, this equation is regularized with a dissipative term, such as viscosity 
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or hyper- viscosity. Additionally Burgers equations can be regularized with linear 
dispersion resulting in the KdV equations. Recently, work has been done by 
our group and others, investigating a regularization of the Burgers equations 
where the nonlinear term is manipulated by averaging the convective velocity 



where g is an averaging kernel with emphasis on the Hclmholtz filter. In this 
paper we extend this technique to the homentropic Euler and Euler equations. 

The technique used on the Burgers equations has been thoroughly investi- 
gated and been established as a valid shock regularization technique. It has 
been shown that solutions to Equations ([ij exist and are unique f?, '3\ for a 
wide variety of filters. This result has also been extended for a multiple dimen- 
sional version of the equations as well [2 . It has been shown that when the 
initial conditions are C^, the solution remains for all time [2 . Furthermore, 
for bell shaped initial conditions it is proven that as the averaging approaches 
zero (a ^ 0), the solutions to Equation (jlj converge to the entropy solution of 
inviscid Burgers equation fS". There is also convincing evidence that this holds 
true for all continuous initial conditions [5 . These analytical results along with 
multiple numerical results regarding shock thickness and energy decay results 
show that this is a valid shock regularization technique. 

The work on the regularization of the Burgers equation is inspired by and 
related to work done on the Lagrangian Averaged Navier-Stokes (LANS-a) equa- 
tions [TllHlElllOliniinillSlITl]. The LANS-a equations also employ an averaged 
velocity in the nonlinear term and have been successful in modeling some tur- 
bulent incompressible flows. 

It is thought that a regularization could be accomplished for the equations 
that describe compressible flow. Encouraged by the results for Burgers equation, 
the next step is to attempt to introduce averaging into the one-dimensional 
homentropic Euler equations, a simplified version of the full Euler equations, 
where pressure is purely a function of density. There have been several attempts 
at such a regularization. 

Inspired by the existence uniqueness proofs from the averaged Burgers equa- 
tions [Ills], Norgard and Mohseni [15 averaged the characteristics of the home- 
ntropic Euler equations to derive the equations 
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with d} — 7p^~^. While these equations were proven to have a convenient exis- 
tence and uniqueness proof, the equations were ultimately found to have signif- 
icant departures in behavior from the homentropic Euler equations, specifically 
for the Riemann problem. What was discovered is that when the characteris- 
tics are averaged, there will be no creation or destruction of characteristics. A 
shock in homentropic Euler equations can produce new characteristics in one of 
the Riemann invariants. The new equations does not capture this behavior and 
results in departures from the desired behavior. 

Using a Lagrangian averaging technique Bhat and Fetecau [16] derived the 
following equations 

Pt + {.pu)x = (3) 

1 9 Vx 

wt + {uw)o, - -{u'^ + a'^uDx = (4) 

Z p 

pw = pv - a^PxU^ (5) 

V ~ U — 0?Uxx- (6) 

While the solutions to the system remained smooth and contained much struc- 
ture it was found that the equations were "not well-suited for the approximation 
of shock solutions of the compressible Euler equations." 

Another attempt by Bhat, Fetecau, and Goodman used a Leray-type aver- 
aging [T7] leading to the equations 

Pt + up^ + pux = (7) 

Ut -I- uux + — = (8) 
P 

u = u- o?Uxx. (9) 

with "p = Kp'^. They then showed that weakly nonlinear geometrical optics 
(WNGO) asymptotic theory predicts the equations will have global smooth 
solutions for 7 = 1 and form shocks in finite time for any other value of 7; 
namely 7 7^ 1. 

Additionally in 2005, H. S. Bhat et. al. (TB] applied the Lagrangian averag- 
ing approach to the full compressible Euler equations. Their approach success- 
fully derived a set of Lagrangian Averaged Euler (LAE-a) equations. However, 
these equations were quite long and complicated, and it seemed that numer- 
ical simulations involving these equations would be impractical for real world 
applications. 

This paper examines another attempt at developing a regularization of the 
homentropic Euler and Euler equations with positive results thus far. Section |2] 
details the motivation behind the technique with section [4] detailing the general 
method. Section[3]specifies the averaging kernels that we consider in the regular- 
izations. Section |6] examines the modified homentropic Euler equations looking 
at conservation of mass and momentum, traveling wave solutions, eigenvalues, 
convergence to weak solutions of the original homentropic Euler equations, and 
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finally examining some numerical results. Section [7] examines the same proper- 
ties, but with the modified Euler equations. All is then followed by concluding 
remarks. 



2 Motivation 

The motivation for our technique stems from a simple concept: that nonlinear 
terms generate high wave modes as continuously as time progresses. Consider 
the mechanics behind shock formation and turbulence. The nonlinear convective 
term u- Vu generates high wave modes, by transferring energy into smaller scales 
as time progresses. This nonlinear term is found in the Burgers equations where 
it causes nonlinear steepening resulting in shocks. It is also found in the Euler 
and Navier-stokes equations where is generates high wave modes by tilting and 
stretching vortices [19] . Thus this nonlinear term is cascading energy down into 
smaller and smaller scales. In the 3D Navier-Stokes and 3D Euler equations, 
the energy cascade has a slope of -| until the Kolmogorov scale, illustrated in 
figure [l^. The Burgers equation has an energy cascade slope of -2 until viscosity 
begins to dominate; seen in figure [T]3 [201 HI] • It is by reducing this cascade 
of energy after some scale that we intend to regularize the Euler equations. In 
the Burgers equation the nonlinear term uu^ was replaced with uu^- A low 
pass filtered velocity will have less energy in its high wave modes after the scale 
— . Thus when inserted into the nonlinear term, the energy cascade will be 
lessened. This modification to the Burgers equations was found to result in a 
regularization of the equations. It is our hypothesis that a similar modification 
to the Euler equations will have similar results. 



3 Averaging kernels 

In our previous work, [21 [5] the low pass filters were assumed to have certain 
properties. Specifically the convolution kernels g were assumed to be normal- 
ized, nonnegative, decreasing, and even. This is summarized in table [T] The 
physical rational behind these properties can be found in our previous work [5]. 
For this paper similar properties are assumed for the filters. However, only sec- 
tions |4] and [s] deal with a general filter. All other sections deal exclusively with 
the Helmholtz filter which possesses all these characteristics. The Helmholtz 
filter is defined as 

u = u-a'^Uxx-, (10) 
and thus has an averaging kernel of 

5"(x) = i-exp(-M). (11) 
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Figure 1: Schematics of energy cascade for the Navier-Stokes/Euler equations 
and the viscous/inviscid Burgers equation are shown. In an inviscid flow, both 
turbulence and shocks show continuous generation of higher wave modes in- 
definitely, (a) Energy cascades from high wavelengths to lower wavelengths at 
a predicted rate of -| for the Navier-Stokes/Euler equations. For the Navier- 
Stokes equations, the kinetic energy drops drastically upon reaching a certain 
wavelength; the Kolmogorov scale, f]. For the Euler equations the cascade con- 
tinues on indefinitely. Here Rei < Re2 < Re-^. (b) For the Burgers equation, 
shocks can form from continuous initial conditions. The energy cascade has a 
slope of -2 when shocks form, until viscosity begins to exert its influence and 
balance the steepening effect of the nonlinear term. Here v^i < U2 < vi- 



Table 1: This table succinctly lists the requirements of the low pass filters 
employed. 



Properties 


Mathematical Expression 


Normalized 
Nonnegative 
Decreasing 
Symmetric 


/.9=1 
5(x) > 0, Vx 
|xi| > |x2| ^ .g(xi) < g(x2) 
|xi| = |x2| =^ .g(xi) ,g(x2) 



4 A general method 

The successful regularization of the Burgers equations can be extended in sev- 
eral ways. As mentioned in the introduction, we had previously interpreted the 
CFB equations as an averaging of the characteristics [TS]. Bhat, Fetecau, and 
Goodman interpreted the regularization of the Burgers equation as an averaging 
of the convective velocity, a Leray-type averaging [17] . Both of these interpre- 
tations, when extended to the homentropic Euler equations, have not led to 
the desired regularization of homentropic Euler equations. Here we examine 
yet another interpretation of the CFB equations, based on a conservation law 
perspective. This approach addresses the cascade of energy generated by the 
nonlinear terms. We then discuss how the technique used in regularizing the 
Burgers equation can be extend to a general technique to be used on conserva- 
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tion laws. Begin by looking at the inviscid Burgers equation 

ut + (uu)^ = 0. (12) 

The flux term represents that the quantity u is flowing into a control volume 
with velocity u. When the product rule is applied to the flux term the result is 

Ut + UUx + UUx — 0. (13) 

Each nonlinear term in the above equation results in steepening waves which 
cascades energy into higher wave modes. It is this energy cascade that produces 
shocks. In section [2] it was discussed how this cascade of energy can be reduced 
by filtering the convective velocity. Thus the non-differentiated term is passed 
through a low pass filter. The resulting equation is 

Ut + UUx + UUx = 0, (14) 

which has been referred to as the convectively filtered Burgers (CFB) equations. 
The portrayal of the CFB equations here differs from Equation ([T]) and in pre- 
vious papers [Il[2lllllll|51IS] by a factor of 2, but is identical under rescaling. 

Now we extend this technique to a more general method. Suppose there is 
a single or multiple conservation laws of the form 

Qt + {Qu)x = 0. (15) 

The proposed method is to apply the product rule to the nonlinear term and 
then apply a filter to the non-differentiated quantities. This results in 

Qt + Qux + uQx = 0. (16) 

We believe that this filtering of the non-differentiated quantities will reduce the 
energy cascade and provide a regularization of the conservation law. It is this 
general method that we apply to the homentropic Euler and Euler equations. 



5 Conservation of conserved quantities 

Before we begin the examination of the modified homentropic Euler and Euler 
equations we will first examine how the modification of the nonlinear term still 
preserves the original conserved quantities. Again consider the more general 
conservation law 

Qt + {Qu)x = (17) 



and its regularized version 



Qt + Qux + uQx = 

u = g * u 

Q^ 9*Q- 



(18a) 
(18b) 
(18c) 
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Integrate Qt over the spatial domain and substitute Equations (18b) and (18c 
and the full definition of convolution to obtain 

^ J Q{x)dx = - j J g{x-y)Q{y)u'{x)dydx- J J g{x - y)u{y)Q' (x) dy dx. 

(19) 

Integrate by parts to obtain 



I J Q{x)dx = J Jg' 



{x - y) {QivHx) + uiy)Q{x)) dy dx. (20) 



If g is even, then g' is odd and g'{x — y) is anti-symmetric over y = x. Clearly 
{Q{y)u{x) + u{y)Q{x)) is symmetric over y = x and thus by symmetry 

^ ^Q{x)dx^Q. (21) 



dt 

From this it is determined the modified conservation law still conserves the 
original conserved quantity. Thus we see that the proposed method conserves 
the quantities that the original conservation laws were designed to preserve. 
This property is independent of the filter so long that it is even, which was a 
requirement of section [3] 

6 The one-dimensional homentropic Euler equa- 
tions 

Examine the homentropic Euler equations 

Pt + (pu)^ = (22a) 
{pu)t + {puu + P)^ = 0, (22b) 



with pressure defined as P = p'' . Equation (22a I and (22b) are the mathe- 
matical expressions of conservation of mass and momentum respectively. The 
homentropic Euler equations makes the assumption that the entropy is constant 
throughout the entire domain and thus the pressure can be expressed solely as 
a function of the density This differs from the isentropic Euler equations 
where the entropy is constant along streamlines but not necessarily constant 
on the whole domain. It is found that the homentropic Euler equations are an 
accurate predictor of gas dynamics behavior for low pressure. 

Both equations have a nonlinear term that fits the general method from 
section [4j Apply the method, break the nonlinear term apart, and then apply 
the filter to the non-differentiated terms in the nonlinear terms to obtain the 
equations 

pt + pu^ + upx = (23a) 
{pu)t + {~pu)ux + u{pu)x + Px^Q (23b) 

with pressure defined as P — p'^ . These equations are now referred to as the 
regularized homentropic Euler equations. For the following analysis, the only 



filter that is considered is the Helmholtz filter (10) 
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6.1 Conservation of mass and momentum 



The homentropic Euler equations ( 22 1 are conservation laws for mass and mo- 
mentum. These equations are reflections of important physical principles. In 
order to be physically relevant it is desirable that any regularization of the ho- 
mentropic Euler equations would also preserve this structure. In section |5] we 
established that the technique used will preserve the conservative structure of 
the original equations. Here we reverify that result for the regularized home- 



ntropic Euler equations ( 23 1 by casting the equations into a conservative form 



using the Helmholtz filter. 

Consider the regularized homentropic Euler equations 



Pt + pux + upx 
{pu)ux + u{pu)x + Px 



with 



P = P 



pu = pu-a pu^^. 



(24a) 
(24b) 

(25a) 
(25b) 
(25c) 



By substituting in Equations ( 25a I , ( 25b I , and ( 25c I and then regrouping 



the similar terms the regularized homentropic Euler equations can be rewritten 
as 



Pt + (pu - a^{upxx + puxx) + a^UxPx).^ = 



{pu)t + (^puu - a^{u{pu),^,^ + puuxx) + a^Uxpu^ + 



P\ =0. 



(26a) 
(26b) 



This shows that the regularized homentropic Euler equations can be written 
in a conservative form and preserve both mass J p and momentum J pu as 
well as the averaged mass / p and averaged momentum J /nZ. Examination of 
Equations (26) could lead to more geometric structure. They may also lead 



to the application of numerical techniques designed specifically for conservation 
laws. 



6.2 Traveling wave solution 

In this section, we establish that a previously known traveling wave solution 
to the homentropic Euler equations is in fact a traveling wave solution to the 
regularized homentropic Euler equations as well. The traveling wave solution 
that we examine for the homentropic Euler equations is a single traveling shock. 

First we establish notation used in this section. The operator [•] will be used 
to quantify the difference between the right and left limits of a function at a 
discontinuity. For example if the limits of u{x) at x* are defined as 

lim u{x) = uji (27) 

X — + 

lim u(x) — ul, (28) 

X — 'X*^ 
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then at X* 

[u] = ur-ul- (29) 

The Rankine-Hugoniot jump conditions establish that for weak solutions of 
conservation laws of the form 

Qt + (f)x = (30) 
a traveling discontinuity must have a speed 5", where S is defined as [22] 



5* = 



[Q]' 



(31) 



Using this relationship it is easy to establish that the following is a weak solution 
to the homentropic Euler equations ( 22 I 

(32) 
(33) 

The speed of the shock S is twice defined, once for each conservation law, by 
the Rankine-Hugoniot jump conditions as 



UL 


X 


< St 


ur 


X 


> St 


PL 


X 


< St 


PR 


X 


> St 



S- 



and 



[p] 

[puu + P] 
[pu] 



Thus ul,ur, and pR must satisfy the relationship 

[pu] ^ ^ ^ [puu + P] 



[P] 



[pu] 



(34) 



(35) 



(36) 



Next we will establish that Equation (32 1 is also a weak solution to the 
regularized homentropic Euler equations ( 26 ) . It is straightforward to establish 
that with u defined as in Equation ( 32 1 and using the Helmholtz kernel definition 
Equation (111 that u and p are 



P = 



exp(^ 


-^*) + 

a > 


UL 


X 


< St 


"^-"« exp( 


x-St\ 


+ Ur 


X 


> St 


P'^-P'^ exp(^ 


-^*)-f 


PL 


X 


< St 


'"'-P"- exp( 


x-St\ 
a f 


+ PR. 


X 


> St 



Thus the values of u and p a,t x = St are 

Ur + Ul 



u{St,t) 
p{St,t) - 



2 

PR + PL 



(37) 
(38) 

(39) 
(40) 
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Additionally using the definition of the Helmholtz filter, u — u — a^Uxx, we 
find that 

(41) 
(42) 



[u\ ^ [u\- a [Uxx] 
Since u is continuous we establish that 

. , 1 . , 



with the similar result 



Ux 



[Px 



u 



(43) 



Looking at the conservation form of the regularized homentropic Euler equa- 
tions ( 26 1 , the Rankine-Hugoniot jump conditions state that in order for Equa- 
tion (|32| to be a weak solution, that the speed of the discontinuity, S, must 

[pu - a^{upxx + puxx) + a^UxPx] 



satisfy 



S ■ 



\P\ 



and 



[puu - a'^{u{pu)^^ + puuxx) + a^Uxpu^ + P] 
[pu] 



Since u,p,Ux, and px are continuous Equation (44 1 reduces to 

-a'^{u[pxx] + p[uxx]) 



S- 



[P] 



Substitute in Equations ( 42 ) and p3| to simplify further to 

iu[p]+p[u]) 



S = 



[P] 



(44) 
(45) 

(46) 
(47) 



^{uji + ul){pr - pl) + Upu + pl){ur - ul) 



[P] 



Now, by using Equations (39) and ( |40| and the definition of [•] one can obtain 

(48) 
(49) 
(50) 



{pbub^ - Plul) 



[P] 



[p] ' 



which is identical to Equation ( 34 1 . Similarly Equation ( 45 ) reduces to Equation 
(35). Thus for this specific example, the Rankine-Hugoniot jump conditions axe 
identical for the homentropic Euler equations and the regularized homentropic 
Euler equations. This validates our claim that Equation (32) is a traveling weak 
solution for the regularized homentropic Euler equations. 
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6.3 Shock Thickness 



With Equation ( 32 ) validated as a traveling weak solution for the regularized ho- 
mentropic Euler equations, we establish an analytical result about shock thick- 
ness. Here the thickness of the shock is defined to be the length over which 90% 
of the amplitude change takes place, centered at the point of inflection. 



Again we note that from Equation ( 32 1 , and the kernel definition of the 



Helmholtz filter (111, the traveling wave solution is 



P = 



"«2"^ exp(^ 


OL / 


UL 


X 


< St 


exp( 


x-St\ 
a z 


+ Ur 


X 


> St 


f^-P"- exp(^ 


a ' 


PL 


X 


< St 


''^7" exp( 


x~St\ 
a f 


+ PR 


X 


> St 



(51) 
(52) 



The thickness of the shock will then be 2ab, where b is the value where 



1 



exp{—\x\) dx 



1 -\x\ 

— exp{^)dx^Q.Q. 
-ab 2a a 



(53) 



This length is independent of pr,pl,ur and ul- As such, the thickness of the 
shock varies linearly on the parameter a. 



6.4 Diagonahzation and eigenvalues 

The homentropic Euler equations have very clearly defined eigenvalues, w ± a, 
where — 'yp^~^ is the speed of sound. This section examines how the pro- 
posed averaging affects the eigenvalues of the system and thus the characteristic 
speeds. In order to cast the equations in vector matrix form we first rewrite the 
equations in their primitive variable form 



pt + pux + upx = 

pu ~ up \ IP''^^ 
Ut + UUx + I I Ux H — Px = 0. 



(54a) 
(54b) 



The equations can then be written in vector matrix form leading to 



P 
u 



— u - 
p 





p 




u 





= 0. 



(55) 



The eigenvalues of matrix A are 




(56) 
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Examining quantities f3 and ^ it seems apparent that as the filtering decreases 

/3 — > and ^ —f 1, thus regaining the original eigenvalues. 

The matrix A can as be diagonalized and can thus be written in the form 



(57) 



where A is a diagonal matrix with its diagonal entries being A^. Using this 
we can write the regularized homentropic Euler equations in the characteristic 
form 

dt dx 



where 



dp 
du 



dv = 

By diagonalizing A we find that 



(58) 
(59) 



2 a? 



pP- ^pjpP'^ + ^a^ p) 
2a2 



(60) 



And thus we can define the quantities through the relation 
dv^ 



2a2 



(61) 



We can then say that dv^ = along the characteristic dx = A* dt. 

If a = 0, there would be no filtering and /3 = 0. In this case we would get 



dv — dp zL —du ^ along the characteristic dx ~ uzLadt. 
a 



(62) 



which is the case for the homentropic Euler equations. For the homentropic 
Euler equations the characteristic variables are able to be computed analytically 
and are 



V — u± 



2a 
7~- 



(63) 



For non-zero a's there appears to be no straightforward way of integrat- 
ing Equation (61 1. Thus currently we have no analytical expression for the 
characteristic variables, though we have not spent much time investigation this 
possibility. 



6.5 Convergence to a weak solution 

It is the goal of our technique to develop new equations that effectively capture 
the low wave mode behavior of the original equations. Thus it is desirable that 
the regularized homentropic Euler equations approximate the homentropic Eu- 
ler equations well. Ideally we can show that as the amount of filtering decreases 
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we regain the original equations. One crucial step in determining this is proving 
that as a —^ the solutions to the regularized homentropic Euler equations 
converge to a weak solution of the homentropic Euler equations. This, among 
other things, proves that as a — > 0, the shocks produced by the regularized ho- 
mentropic Euler equations travel at the same speed as those of the homentropic 
Euler equations, which is desirable. 

As of yet, we have not been able to extend our existence theorems of the reg- 
ularized Burgers equations ([T]) to the regularized homentropic Euler equations 
(23). Without this theorem we are forced to make several assumptions that 



we consider modest, considering the numerical evidence presented in section 



6.7 We assume that for every a > there exists a solution to Equations (24 1 



Beyond that it is assumed that a subsequence of those solutions converge in 
^loc ^^'^ solutions are bounded independent of a. The following summarizes 
these assumptions. 



|w||oo < U 

IIpIIoo <i? 

„ = u in Ll^ 
lim p = pin 

a— *0 



lim u ~ 



loc 



(64a) 
(64b) 
(64c) 

(64d) 



With these assumptions we are able to prove that the solutions to the regu- 



larized homentropic Euler equations ( 23 1 will converge to weak solutions of the 



homentropic Euler equations ( 22 1 . The examination of the claim is done with 
the Helmholtz filter, with u defined as 



with 



and thus 



u = u * g 



g{x) = -exp{-\x\) 



g«(x) = ;^exp(-M). 
2a a 



(65) 



The unfiltered and filtered velocities can also be related by 



u = u — a Ur 



(66) 



Both definitions are used. The following bounds are easily established by ex- 
amining the kernel. 

Ilfflli = l 

and 

9 =1, 



and thus 



dx 
d 

Tx' 



1 

a 
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These bounds combined with Young's inequaHty "23] and the assumptions ( 64 1 
give the following estimates 



||u||oo < U 

IIpIIoo < R 

\pu\\oo < UR 



< -R 



< -UR 



(67) 
(68) 
(69) 

(70) 
(71) 
(72) 



Along with these estimates, the last piece needed is taken from Duoandikoetxea 
[24l . The following lemma is a restatement of Duoandikoetxea's Theorem 2.1 
from page 25. 

Lemma 6.1 Let g be an integrahle function on M such that J g = 1. Define 
9"=h9iV- Then 

lim||5"*/-/||p = 

a— >0 

if f Cz < p < oo and uniformly (i.e. when p — oo) if f £ Co(II^)- 

With lemma [6T| the convergence to weak solutions can now be proven. Begin 



by multiplying Equations ( 26 1 by a test function <j) and integrate over time and 
space. It is assumed that has an infinite number of bounded and continuous 
derivatives and is compactly supported. Doing this we obtain the equations 



[pu - a^{up^x + puxx) + a^UxPx)^4'dtdx = 

(73a) 

{pu)t(t>+ (jMu - a^{u{pu),^^ + puuxx) + a'^UxpUx^ + P)x4>dtdx = 0. 

(73b) 



Integrate by parts to obtain 



bt + {pu) (j)x dtdx 







{pu)(j)t + {puu + P) 4>x dt dx 



{a^{upxx + puxx) - a^UxPx) 4>x dt 
(74a) 

(^a^{u{pu)^^ + puuxx) - a'^UxpUx))^xdtc 

(74b) 



Clearly if the right hand side of Equations ( 74 1) and ( 74 d) limit to zero then the 
limit of p and m is a weak solution to the homentropic Euler equations. Begin 
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by examining the term 



a^upxx(j)x dt dx 



(75) 



Substitute the definition of the averaging so that o? Pxx = P ~ P and examine 
the quantity 



u{p — p)4>x dt dx. 



(76) 



It was estabhshed that u and 4>x are bounded and by lemma 6.1 {p— p) converges 
in to zero. Thus we find 



hm 



a upxx<t>x dtdx ■ 



u{p ~ p)(l)x dt dx — 0, (77) 



and the similar terms can be treated likewise. 
Next examine the term 



a Uxpu^(j)x dt dx. 







Perform an integration by parts to obtain 

r-T 

a'^u{pu)^^(j)x dtdx - 



a upu^4)xx dt dx. 



(78) 



(79) 



The first term in Equation (79 1 limits to zero from the steps shown above. The 
second term can be bounded with the estimates found from Young's inequality. 



a^upu^4>xx dt dx 







< 



a^U-UR\(t>xx\dtdx (80) 
a 



^aU^RW^xxWi (81) 

and thus limits to zero. The similar terms can be treated likewise. Thus we see 
that solutions to Equations ( 26 1 converge to weak solutions of the homentropic 
Euler equations as a — > 0. 



6.6 Numerics 

In this section, we discuss the numerical techniques we used to simulate the 
regularized homentropic Euler equations and the homentropic Euler equations. 
Since the regularized homentropic Euler equations appear to be regularized and 
its solutions are smooth, we are able to use a pseudo-spectral method. The 
homentropic Euler equations tend to form shocks, so we use well-established 
techniques to handle this behavior. 
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6.6.1 Numerical simulations of the regularized homentropic Euler 
equations 

For the regularized homentropic Euler equations, much like in turbulence simu- 
lation, we wish only to resolved the averaged quantities in our numerical simu- 
lations. To obtain the equations with only the smooth variables, we first apply 



the Helmholtz filter to Equations (23 1 and add {pu)x or {puu)^ to both sides. 



Pt + {pu)^ = {pu)x - {pux + upx) (82a) 
(OUj -I- (puw + P)^ = iW^)x ^ iWUx + u{pu)x). (82b) 

Then use the definition of the Helmholtz filter ( |To| ) to manipulate the equations 
to 



Pt + (pu)^ = (1 - a^dl) [pu)^ - p{ux - a^Uxxx) - u{Px - a^Pxxx) 

(83a) 

pUf + {puu + P)^ = (1 - a^dl) {puu)x - Wiux - a^Uxxx) - u{{pu)x - a'^{pu)xxx)- 

(83b) 

One can then simplify to obtain 



Pt + {pu)x = -3q! {uxPx)x (84a) 
puj + {puu + P)^ = ~'ia^{uxpux)x- (84b) 

The equations are now fully in terms of the averaged quantities with the right 
hand side being the regularizing terms. 

In order to reduce the potential for numerical dissipation we use a pseudo- 



spectral method to solve Equations ( 84 1 . We advance the equations in time with 



an explicit, Runge-Kutta-Fehlberg predictor/corrector (RK45) [6^. The initial 
time step is chosen small enough to achieve stability, and is then varied by the 
code using the formula 

h^+.^lhJ^-^Krr-X ■ (85) 
\\\Pt - pAhJ 

Thus the new time step is chosen from the previous time step and the amount of 
error between the predicted density, p and the corrected density p. The relative 
error tolerance was chosen at e = 10""* and the safety factor 7 = 0.9. If the 
new time step chosen was found to violate the CFL condition the time step 
was chosen according to the velocity speed and the speed of sound with CFL 
number 0.5. 

Spatial derivatives and the inversion of the Helmholtz operator were com- 
puted in the Fourier domain. The terms were converted into the Fourier domain 
using a Fast Fourier Transform, multiplied by the appropriate term and then 
converted back into the physical domain. 

A number of high resolution simulations were done at the resolution of 
2^"* = 65536 grid points on the domain [0,27r]. We found that some long term 
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simulations may suffer from long term numerical instabilities. After approxi- 
mately 10,000 time steps, we would notice energy building up in the high wave 
modes near the resolution boundary. If not checked, these would contaminate 
the results of lower wave modes. In order to control this long term instability, 
every 200 time steps we zero out all the wave modes with wave number higher 
that y where N is the highest Fourier wave mode simulated. Wave modes 
higher than are set to zero at every nonlinear multiplication to prevent 
aliasing. Numerical runs were conducted for values a = 0.10, 0.09, 0.02, 0.01. 
In the worst case scenario, a = 0.01, the wave modes that are zeroed are over 
an order of magnitude higher than ^ , as seen in figure |2] 



N/3 2N/3 




Figure 2: The spectral energy of u. The vertical lines represent three significant 
wave numbers. From left to right the lines represent y, and where N 
is the highest Fourier wave mode simulated. The wave mode y is the highest 
wave mode that is ever zeroed out in the simulations. Here with a = 0.01 you 
can see that zeroing of wave modes takes place an order of magnitude higher 
than where the filtering has its effect. This graph was taken from simulations 
of examples problem 90 Here t = 0.1, N = 2^^, and n = 2^"*, where N is the 



number of Fourier modes resolved and n is the number of grid points. 



We attribute this error to the right hand side of Equations (84) for the 
following reasons. Examine the term 



- 3a {u^,px)x- (86) 

Expand the term to 

3q; UxxPx ^xPxx- 

(87) 

and specifically notice the UxPxx term. This term can be considered as a 
viscosity- like term that changes signs with Ux- Now along the expansion wave 
Ux > Q and thus there is a negative second derivative or viscosity-like term. 
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Negative viscosity is inherently unstable, so we attribute the numerical insta- 
bilities to this. This is similar to the backscatter of energy to larger scales that 
is observed in turbulence. 



6.6.2 Numerics for the homentropic Euler equations 

For the numerical simulations of the homentropic Euler equations there are 
many established, available techniques |22j. We chose to use the Richtmyer 
method, a well-established if low order method, was utilized as described by 
|22l. This method is a second-order, finite-difference scheme and employs an ar- 
tificial viscosity. This method requires an artificial viscosity for stability when 
examining the Riemann problem. Several different values of v were tested to 
see that the value did not significantly affect the solutions on the time interval 
examined. For the numerical simulations shown here, the artificial viscosity was 
set at J/ = 0.08. For reference, the simulations were done with 2^'* grid points 
on a [0, 27r] domain. 



6.7 Numerical results 



This section examines some numerical simulations performed on Equations ( 84 1 
with the technique described in the previous section. A shock tube or Riemann 
problem exhibits both shocks and expansion waves, some of the key behaviors of 
the homentropic Euler equations. Our pseudo-spectral method enforces periodic 
boundary conditions, so we created two pressure jumps in the initial conditions 
resolving this issue. Thus essentially a double shock tube problem is considered 
where there are two discontinuities in the initial conditions to make the left and 
right side boundary conditions identical. The example problem considered is 

uq(x) = (88) 

r 1 o<x<f 

ao{x) = I 2 ^ <x<f (89) 
[ 1 <x<2n 

which for 7 = 1.4, the constant for air, is equivalent to 

uo{x) = (90) 

r 0.43120 <x<^ 
po{x) = I 13.7984 ^ <x<^ . (91) 
[ 0.43120 ^ <x<2t: 

In our previous work, we established that for the CFB equations ([T|) initial 
conditions with discontinuities were excluded in order to avoid non-cntropic be- 
havior [5j . By regularizing the initial conditions it was proven that the solution 
to the CFB equations would be regularized for all time. This was done by av- 
eraging the initial conditions with the same filter used on the velocity. We use 
this same approach here. Thus when the filter is applied to the entire equations 
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to obtain Equations ( 84 ) , the initial conditions arc filtered twice with the same 



filter used on the equations. 

In the simulation of the double shock tube problem the regularized homen- 
tropic Euler equations are displaying similar behavior to the homentropic Euler 
equations. Figures [3] and |4] show the solutions to the two sets of equations 
imposed on each other. One can see that the solutions to the regularized home- 
ntropic Euler equations capture both the expansion wave and the shock front. 
The solutions to the regularized homentropic Euler equations are seen to be 
smooth with the solutions tightening to the solutions of the homentropic Euler 
equations as a decreases. 





Figure 3: This figure shows a numerical simulation of the regularized homen- 
tropic Euler equations (dashed line) plotted against the solution to the homen- 
tropic Euler equations (solid line). Here the value of a = 0.05. In both figures, 
it is clear that the regularized homentropic Euler equations are capturing both 
the expansion wave and shock behavior, (a) The velocity, (b) The density. 



Now we check the convergence of the solutions of the regularized homentropic 
Euler equations to the solution of the homentropic Euler equations as a ^ 0. 
Figures [s] and |6] show that as a ^ the error in the norm appears to be 
approaching zero for the example problem. This suggests that the solutions of 
the regularized homentropic Euler equations converge to the solutions of the 
homentropic Euler equations. 

6.8 Kinetic energy rates 

In addition to checking solution profiles, we examine the effect that our reg- 
ularization technique has upon the kinetic energy. For the homentropic Euler 
equations we define kinetic energy as | pu^ . For the regularized homentropic Eu- 
ler equations there are three different averaged quantities, p, u, and pu. With 
these quantities, kinetic energy can be defined in a variety of ways. In this 
section we examine a kinetic energy with unfiltered terms, -^pu^, and a kinetic 
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Figure 4: This figure shows a numerical simulation of the regularized homen- 
tropic Euler equations (dashed line) plotted against the solution to the homen- 
tropic Euler equations (solid line). Here the value of a = 0.01. In both figures, 
it is clear that the regularized homentropic Euler equations are capturing both 
the expansion wave and shock behavior. With the lower value of a the fit is 
much closer, (a) The velocity, (b) The density. 




Figure 5: This figure displays the difference between the density in solutions 
of the regularized homentropic Euler equations and the solution of the ho- 
mentropic Euler equations in the norm as a — )■ 0. As a — > the differ- 
ence in the solutions also approaches zero. The measurements were taken for 
a = 0.01, 0.02, 0.1 at times t = 0, , t = 0.2 , and t = 0.4 . 



energy with filtered terms, \p'>j? . We have also examined the possible kinetic 
energies -^^u and but for the example we are considering the difference 
between them and \pv? was of the order 10~^^ and thus negligible. For flows 
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0.02 0.04 0.06 0.08 o!l 

a 

Figure 6: This figure displays the difference between the momentum in solu- 
tions of the regularized homentropic Euler equations and the solution of the 
homentropic Euler equations in the norm as a — *■ 0. As a — *■ the differ- 
ence in the solutions also approaches zero. The measurements were talcen for 
a = 0.01, 0.02,..., 0.1 at times t = 0.2 , and t = 0.4 . 



with more small scale behavior, we would expect this difference to be more 
significant. 

For the shoclc tube problem, the fcinetic energy of the system clearly starts at 
zero. For the homentropic Euler equations the solution is self similar, depending 
only on the variable j. Thus the kinetic energy for the homentropic Euler 
equations will be a linear function of time. This can be seen in figure [7] We 
find that the energies of the regularized homentropic Euler equations mimic 
that of the homentropic Euler equations. In the simulations of the regularized 
homentropic Euler equations there is a brief period where the energy growth is 
curved before it appears to behave linearly. We attribute this to the averaging 
of the initial conditions. As we would expect, as a decreases the energies grow 
closer to those of the homentropic Euler equations. 

7 The one-dimensional Euler equations 

While the homentropic Euler equations are good for low pressures, to have real 
impact we want our technique to be able to capture the behavior of the full one- 
dimensional Euler equations. With the regularized homentropic Euler equations 
showing promise, we attempt to use the same general methodology on the one- 
dimensional Euler equations. The one-dimensional Euler equations consist of 
three conservation laws paired with a constitutive law. The conservation laws 
are conservation of mass, momentum, and energy. For the constitutive law, this 
paper is considering an ideal gas. Thus the one-dimensional Euler equations 



^1 
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Figure 7: The kinetic energy of the homentropic Euler equations and the reg- 
ularized homentropic Euler equations. The energy for the true solution of the 
homentropic Euler equations is shown as a dashed red line. The simulation of 
the regularized homentropic Euler equations for a = 0.01, 0.02, ...,0.1 are shown 
as solid black lines. The bottommost line represents a = 0.1. As a decreases the 
energy approaches the energy of the homentropic Euler equations, (a) These 
are plots of the unfiltered kinetic energies, ^pu^. (b) These are plots of the 

filtered kinetic energies, When examined the plots of ■^'puu and were 

identical to this one. 



considered, in conservation form, are 



Pt + (pu)^ = (92a) 

{pu)t + {puu + P), - (92b) 

{pe)t + {peu + uP)x ^ (92c) 

P={j-l)(pe^l-puA (92d) 



The technique described in section |4] applies the product rule to nonlinear 
terms and then the non-differentiated quantities are spatially filtered. Again this 
is to reduce the production of higher wave modes as time progresses, although 
as in the regularized homentropic Euler equations, we will see that an artificial 
viscosity is needed for numerical stability. When our averaging technique is 
applied, the new equations are 

Pt + pux + upx = (93a) 
{pu)t + puux + u{pu)x + Px ^0 (93b) 
{pe)t + peux + u{pe)x + Pu^ + uPx = (93c) 

P = (7 - 1) (pe - ^puA (93d) 
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which we refer to as regularized Euler equations. The general method estab- 
lished in section |4] does not address exactly how to handle the pressure terms. 
In the conservation of momentum equation ( 92b ) the pressure term is left unaf- 
fected as in the regularized homentropic Euler equations. However, the {uP)^ 
term in the conservation of energy equation ( 92c I is averaged using the general 
method. We have found in earlier numerical simulations that not performing 
the averaging technique on this term lead to some nonphysical behavior. 

The rest of this paper is dedicated to the examination of the regularized 
Euler equations. The next section shows that the regularized Euler equations 
conserve mass, momentum, and energy, by casting the equations into a con- 
servative form using the Helmholtz filter. Section |7.2| establishes a traveling 



wave solution using the same techniques as before. Section 7.4 proves that with 
modest assumptions the regularized Euler equations converge to a weak solution 
of the Euler equations as a — > 0. Finally sections |7.5| and |7.6| examine some 
numerical simulations and their results. 



7.1 Conservation of mass and momentum 



The Euler equations (92) are conservation laws that address conservation of 
mass, momentum, and energy. This section show that the regularized Euler 
equations preserve these same quantities by casting the equations into a con- 
servative form. It would be sufficient to refer the reader back to sections [5] and 
16.11 and note that our method does not disturb the conservation structure of 
the Euler equations. However, explicitly showing this for the regularized Euler 
equations for the Helmholtz filter has merit as it demonstrates a new form of 
the equations which arc of interest. 

To review, the Helmholtz filter the non-filtered quantities can be expressed 



by their filtered counterparts by Equation ( 10 1 . Using this filter. Equations ( 93 1 
can be rewritten as 



pt+ [pu - a'^iupxx + puxx) + a^UxPx]^ = 



{pu)t+ puu - a'^{u{pu)^^ + puuxx) + a^Uxpu.^ + P 
{pe)t+ peu - a^{u{pe)^,^ + peuxx) + a'^Uxpe^+ 
+ Pu- a^{uP^^ + Puxx) + a^UxPx] ^ = 



(94a) 
(94b) 

(94c) 



In this conservation form it is easy to see that mass J p, momentum J pu, and 
energy J pe are conserved. 



7.2 Traveling wave solution 

In section |6.2| it was established that certain traveling weak solutions to the 
homentropic Euler equations were also solutions to the regularized homentropic 
Euler equations. The section establishes a similar result for the Euler and 
regularized Euler equations using the same techniques. Again the operator [•] 
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will be used to quantify the difference between the right and left hmits of a 
function at a discontinuity so that 



[u\ =ur- ul- 
Again we examine a traveling shock of the form 

P = 
pu ~ 
pe = 



PL 


X < St 


PR 


x> St 


Ml 


X < St 


Mr 


x> St 


El 


X < St 


Eb. 


x> St 



(95) 

(96) 
(97) 
(98) 



Using the Rankine-Hugoniot jump conditions (31 1 we can establish that at a 



discontinuity a weak solution to the Euler equations ( 92 1 must satisfy 

[pu] 



S = 
S = 
S = 



[P] 

[puu + P] 

[pu] 
[peu + uP] 
[pe] 



(99a) 
(99b) 
(99c) 



The Rankine-Hugoniot jump conditions for the regularized Euler equations 
are 



S = 
S = 
s = 



[pu - a'^{upxx + puxx) + a^UxPx 
[P] 



(100a) 



(100b) 



\puu - a^{u{pu)^^ + puu,j,^) + cPu^pu^ + P] 
[pu] 

[peu ~ a'^ju^pe)^^ + peUxx) + a^Uxpe^ + Pu - a'^{uPxx + Puxx) + ^^Yf f^A 



Using the exact same computations as in section 6.2 it is a straight for- 
ward process to show that for solutions of the form (96 1, the jump conditions 



for the Euler equations (99 1 are exactly equivalent to the jump conditions for 



the regularized Euler equations (100 1. Thus Equation (96) with values satis- 
fying Equations ( 99 1 will be traveling weak solutions to the regularized Euler 
equations. 

Additionally using the same analysis found in section |6.3| it can again be 
showing that the thickness of the shocks for the traveling solution will decreases 
linearly with a. 
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7.3 Eigenvalues 



Much like the homentropic Euler equations, the one-dimensional Euler equations 
have very well defined eigenvalues u, u±a. Again we examine how the averaging 
technique has affected the eigenvalues of the system. In this section we examine 
the eigenvalues of the regularized Euler equations. We begin with Equations (93 1 
through much substitution and manipulation you can express them in primitive 
variable form as such: 



Pt 



upx 



pu-x = 
f~pu~ up 



p 



(101) 
(102) 



Pt - (7 - l)PpuUx + (7 



l)P+(7-l)pe--(7-l)P"' 



+ (7U - (7 - l)u)Px = 0. 



The equations can then be written in vector matrix form leading to 



p 




u 


u 


+ 





p 


t 


_ 



p 

13 + u 
-(7 - l)l3pu + pa^ 



7U — (7 — l)u 







p 








u 


= 


l)u _ 




p 


X 



(103) 



(104) 



Define the quantities [3, k, and 5, as 
pu — up 



P 

- u 



((7-l)P+(7-l)pe- 
P 



\{l~l)-pu') 



Then the eigenvalues of matrix A are 
Ao = -0 



(105) 
(106) 
(107) 

(108) 



A+ = 



7K 



± 



4 



(7-1)2^2 (7-l)(u + u)/3 



The quantities (3 and k would appear to limit to zero as a — > 0, and the 
quantity a? would appear to limit to a? , the speed of sound for the Euler equa- 
tions. Again the averaging technique alters the eigenvalues, but with the original 
values regained with the limit as a — > 0. 

Much like in the homentropic case the matrix A can as be diagonalized and 
can thus be written in the form A = QAQ~^. Using this we can write the 
regularized Euler equations in the characteristic form 

dv ^dv ^ 

dt dx 



(110) 
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where 



dp 
du 
dP 



When diagonahzing A we calculated 



a -(7-1)"^ 
(7-1)kP 






P(7-1)k 

(7-l)K-;3+^/32 + (7-l)3K2-2(7-l)(M+u)/3+4a2 
2p(a2-(7-l)M/3) 

(7-l)K-ff-y/3^ + (7-l)^K^-2(7-l)(fi+u)/3+4a^ 
2p(a2-(7-l)«/3) 



(111) 



(112) 



With Q^^ containing such a complex structure we were not able to find a 
straightforward way of determining an analytical expression for the character- 
istic variables. 



7.4 Convergence to a weak solution 



Section |6.5| showed that with certain assumptions made on the solutions of 
the regularized Euler equations, the solutions converge to weak solutions of 
the homentropic Euler equations as a — > 0. With similar assumptions on the 
regularized Euler equations, we are able to prove that the solutions converge to 
weak solutions of the Euler equations. 

Again the importance of showing that solutions of the regularized Euler 
equations converge to weak solutions of the Euler equations is that it is desirable 
that the regularized Euler equations capture the physical behavior demonstrated 
in the Euler equations. Among these is shock speed. Convergence to a weak 
solution verifies that shocks formed in the NME equations have speeds that 
limit to the proper shock speeds found in the Euler equations. 

With an existence and uniqueness proof for the regularized Euler equations 
not yet developed, we are forced to make several assumptions, that are modest 



considering the numerical evidence presented in section 7.5 Assume that for 
every a > there exists a solution. Beyond that assume a subsequence of those 
solutions converge in Lj^^ and the solutions are bounded independent of a. The 
following summarizes these assumptions. 



lim u 
lim n 
lim e 



\\u\\oc < U 


(113a) 


\\p\\oo<R 


(113b) 


||e||oo < E 


(113c) 


= u in L}^^ 


(113d) 


= P in Lj^c 


(113e) 


= e in L}^^ 


(113f) 




(113g) 
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With these assumptions we are able to prove that the weak solutions to 



the regularized Euler equations ( 93 1 will converge to weak solutions of the Eu- 
ler equations (92 1 as a 0. The examination of this claim is done with the 
Helmholtz filter which has bounds already established in section |6.5[ These 



bounds can be combined with Young's inequality [27" and assumptions (1131 to 



obtain estimates on the filtered quantities. The following estimates essentially 
state that the filtered quantities have the same infinity bound as the unfiltered 
quantities and that the first derivatives of the filtered quantities are of order — . 
Explicitly these estimates are 



||m||oo < u 

\\p\\oo<R 

\pu\\oo < UR 
||e||oo < E 
\\pe\\oo < EE 

\P\\oo <h~l)[ER+^RU^ 



\Px\ 



< -U 



oc < —R 

a 

< -UR 
a 

< -ER 
a 



\P.\\oo < -(7-1) (eR+Iru'' 
a \ 2 



(114) 
(115) 
(116) 
(117) 
(118) 

(119) 

(120) 
(121) 
(122) 
(123) 

(124) 
(125) 



Begin by multiplying Equations ( 94 1 by a test function and integrate over 



time and space. The test function (j) has an infinite number of bounded and 
continuous derivatives and in compactly supported. After the multiplication 
and integration the equations are now 



Pt0+ {pu - a^{upxx + puxx) + a'^UxPx)^4'dtdx = 



(126a) 



{pu)t(i)+ [puu - a^{u{pu)^^ + puuxx) + a^Uxpu^ + P)xcj)dtdx = 

(126b) 



{pe)t4>+ peu - a'^{u{pe)^^ + peuxx) + o?Uxpe^+ 

+ Pu — o?(uPxx + Puxx) + ct'^UxPx] ^(j)dtdx — 0. 



(126c) 
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Integrate by parts to obtain 



p(j>t + {pu) (p^dtdx ^ I I {a^iupxx + puxx) - a^UxPx) <t>x dt dx 



(127a) 

{pu)4>t + {puu + P) (jj^dtdx ^ / (^a^{u{pu)^^ + ~puuxx) - a^u.^Wi^)4>xdtdx 



T I- i-T 



LI 



(127b) 

{pe)4>t + {peu + Pu) c/jx dtdx = / (^a'^{u{pe)^^ + peuxx) - a^Uxpe^y, 



: dt dx 



T 

{a^{uP^^ + Puxx) - ci^UxPx))4'x dtdx. 

(127c) 







Clearly if the right hand side of the above equations limits to as a ^ then 
convergence to a weak solution is proven. Since the calculations are all but 
identical, the reader is referred to section |6.5| to prove that the right hand side 
of the equations do limit to 0. With the right hand side limiting to we find 
that p, u, and e are, in fact, limiting to a weak solution of the Euler equations. 

7.5 Numerics 

The numerical simulations for the regularized Euler equations are the same as 
those discussed in section |6.6| with the addition of the energy equation. The 
equations simulated here are 



Pt + (pu)^ = - 3a {uxPx)x (128a) 
+ {puu + P) = - ?,a^{ux'pUx)x (128b) 



pet + {-peu + uP)^ = - ?,a^{ux {pCx + Px)).^- (128c) 



The same pseudo-spectral method is used to solve Equations (1281. Again 
an adaptive Runge-Kutta is used to advance in time with all spatial derivatives 
and the inversion of the Helmholtz operator conducted in the Fourier domain. 
The simulations were conducted with 2^^ = 65536 grid points on the domain 
[0, 27r] for approximately 10,000 time steps. Numerical runs were conducted for 
values a = 0.10, 0.09, 0.02, 0.01. The same long term instability was found in 
these numerical simulations and were again controlled by setting wave modes 
higher than y to zero every 200 time steps. 

7.6 Numerical results 

Again the preferred example problem is the shock tube or Riemann problem. 
This is chosen for its demonstration of expansion waves, shocks, and contact 
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surfaces. The example problem is one of the classic Sod test problems [5S] 
where the initial conditions are 



uq{x) = (129a) 

^«(-)- { J.125 l<l% (129b) 

o / X r 10 < a; < TT nnn \ 

= [l n<x<2n (129c) 

These initial conditions were chosen as they produce an expansion wave, a con- 
tact surface, and a shock, the three classical behaviors of the Riemann problem 
in the Euler equations. As with the homentropic Euler equations the initial 
conditions are twice filtered. 

The numerical method that we chose forces periodic boundary conditions. 
Thus for the shock tube problem we chose, errors begin to propagate from the 
edges as time progresses. While the simulations were conducted on [0, 27r] we 
only consider the results valid for , ^] up to time t — 0.25. 

Figures |8] and [9] show simulations of the regularized Euler equations for two 
different values of a plotted against the solution to the Euler equations. The 
expansion wave, contact surface, and the shock are all captured in the behavior 
of the regularized Euler equations. For the smaller value of a, the behavior of 
the regularized Euler equations matches the Euler equations more closely. 




234 "'234 ''^234 




Figure 8: This figure shows a numerical simulation of the regularized Euler 
equations (dashed line) plotted against the solution to the Euler equations (solid 
line). Here the value of a = 0.05 at time t — 0.25. In the figures, it is clear that 
the regularized Euler equations are capturing both the expansion wave, contact 
surface, and shock behavior, (a) The density, (b) The velocity, (c) The energy. 



As before we check the convergence of the solutions of the regularized Euler 
equations to the solution of the Euler equations as a ^ 0. Figures [TO] [TT] and 
12 show that as a — > the error in the norm appears to be approaching zero 



for the example problem. This suggests that the solutions of the regularized 
Euler equations converge to the solutions of the Euler equations. 
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Figure 9: This figure shows a numerical simulation of the regularized Euler 
equations (dashed line) plotted against the solution to the Euler equations (solid 
line). Here the value of a = 0.01 at time t — 0.25. In the figures, it is clear that 
the regularized Euler equations are capturing both the expansion wave, contact 
surface, and shock behavior. With the lower value of a the fit is much closer, 
(a) The density, (b) The velocity, (c) The energy. 




Figure 10: This figure displays the difference between the density in solutions 
of the regularized Euler equations and the solution of the Euler equations in the 
norm as a ^ 0. As a ^ the difference in the solutions also approaches 
zero. The measurements were taken for a = 0.01, 0.02, 0.1 at times t = 0.05, 
, t = 0.15 , and i = 0.25 . 



7.7 Kinetic energy rates 

As in section |6.8| we examine the kinetic energy for the shock tube problem. 
Again for the shock tube problem, the solution to the Euler equations are self 
similar, depending only on the variable j and thus the kinetic energy changes 
linearly in time. For the Euler equations, we examine the kinetic energy ^pu^ 
and for the regularized Euler equations we examine an unfiltered kinetic energy 
^pu^ and a filtered kinetic energy ^pu^. Other filtered kinetic energies were 



30 




Figure 11: This figure displays the difference between the momentum in solu- 
tions of the regularized Euler equations and the solution of the Euler equations 
in the norm as a ^ 0. As a — > the difference in the solutions also ap- 
proaches zero. The measurements were taken for a = 0.01, 0.02, 0.1 at times 
t = 0.05, , t = 0.15 , and t = 0.25 . 




Figure 12: This figure displays the difference between the energy in solutions of 
the regularized Euler equations and the solution of the Euler equations in the 
norm as a ^ 0. As a — > the difference in the solutions also approaches 
zero. The measurements were taken for a = 0.01,0.02, ...,0.1 at times t = 0.05, 
, t = 0.15 , and t = 0.25 . 



also considered, but as in the homentropic case, for this example problem, the 
differences between them and -^pu^ were negligible. 

Figure [13] shows how the kinetic energies for the regularized Euler equations 
behave for various values of a. Again after a brief period, the energies seem 
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to vary linearly with time. We attribute this brief period to the averaging of 
the initial conditions. As a decreases the energies for the regularized Euler 
equations approach the energy of the Euler equations, as we would expect if the 
solutions are converging as a — > 0. 




Figure 13: The kinetic energy of the Euler equations and the regularized Euler 
equations. The energy for the true solution of the Euler equations is shown 
as a dashed red line. The simulation of the regularized Euler equations for 
a = 0.01, 0.02, 0.1 are shown as solid black lines. The bottommost line 
represents a — 0.1. As a decreases the energy approaches the energy of the 
Euler equations, (a) These are plots of the unfiltered kinetic energies, ^pu^. 
(b) These are plots of the filtered kinetic energies, ^pu^. When examined the 

plots of ^puu and were identical to this one. 



8 Conclusion 

Using the convectively filtered Burgers (CFB) equations as inspiration we have 
developed a new averaging technique with the intent of regularizing both shocks 
and turbulence simultaneously. This paper examines primarily the shocks regu- 
larization aspect of the technique. We discussed the physical motivation for the 
method and then derived a general technique to be used on conservation laws. 
It was then established that this technique, when applied to conservation laws, 
would preserve the original conservative properties. 

The remainder of the paper then examined the effects when this method was 
applied to the homentropic Euler and Euler equations. The results show much 
promise. It was found that with the Helmholtz filter that both the regularized 
homentropic Euler and regularized Euler equations can be rewritten in conser- 
vation form, reiterating that the original conservative properties are preserved. 
For both sets of equations we were able to find traveling shock solutions, where 
the Rankine-Hugoniot jump conditions for the modified equations reduced to 
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the same jump conditions of the original equations. For both sets of equations 
it was proven that as the filtering is decreased, a ^ 0, the solutions will con- 
verge to weak solutions of the original equations. Both these results show that 
the proper shock speeds of the original equations will be preserved with the 
averaging techniques. 

Numerical simulations were run on both sets of equations for the shock 
tube problem. These simulations demonstrated that the modified equations 
mimic the behavior of the original equations. The regularized homentropic Euler 
equations captured both the expansion wave and the shock front behavior, while 
the regularized Euler equations captured the expansion wave, contact surface, 
and the shock front. The solutions appeared regularized meaning that they were 
continuous and smooth. Furthermore, as a — > these solutions were seen to be 
converging to the solutions of the original equations. 

There is still more work to be done regarding these equations. It would be 
beneficial to establish more theory regarding both sets of equations. Specifically 
existence proofs would be beneficial. It would be interesting to see if the either 
sets of equations possess a Hamiltonian structure. 

The regularized homentropic Euler and regularized Euler equations are show- 
ing promise as a new regularization method based on the preliminary examina- 
tion considered here. From these results we believe that this averaging technique 
leads to a regularization of the homentropic Euler and Euler equations that is 
capable of capturing the relevant behavior of the equations, at least for short 
time simulations. With future work we hope to establish this more thoroughly 
and extend this technique into higher dimensions. 
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